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Abstract 

In this work we investigate the effects of the littlest Higgs model (LHM) up to 
the QCD next-to-leading order (NLO) on the Z°H° associated production at the 
CERN Large Hadron Collider (LHC). We study the dependences of the leading or- 
der and NLO QCD corrected integrated cross sections for this process on the factoriza- 
tion/renormalization scale and the LHM parameters. We also provide the distributions 
of the transverse momenta of final decay products yT and r~ . Our results show that 
the heavy neutral gauge bosons Z# and An could induce significant discrepancies from 
the standard model predictions. It is found that when the LHM parameters are taken 
as c = 0.5, d = 0.22, / = 4 TeV and /j, = (M H +M z )/2, the effects at the y/s = 14 TeV 
LHC from the heavy neutral gauge boson are about 12.83% and 10.37% to the leading 
order and NLO QCD corrected integrated cross sections, respectively. We also con- 
clude that the NLO QCD corrections at the y/s = 14 TeV LHC can obviously reduce 
the scale uncertainty of the integrated cross section, and significantly enhance the dif- 
ferential cross sections of pj, and pTp . It demonstrates that the precision measurement 
of the Z°H° associated production process at the LHC could provide the clue of the 
LHM physics. 



PACS: 12.38.Bx, 14.70.Hp, 14.80.Cp 
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I. Introduction 



The CDF and DO experiments ruled out the standard model (SM) [fl [2] Higgs boson with 
mass between 156 GeV and 177 GeV at 95% confidence level (CL) [3]. Recently, the ATLAS 
and CMS experiments at the CERN Large Hadron Collider (LHC) have excluded most of 
the Higgs mass ranges of 146 — 466 GeV and 145 — 400 GeV at 95% CL in their reports 
of [1] and [5], respectively. Currently, the ATLAS and CMS groups exclude a substantial 
region of the possible Higgs boson mass range, and find several Higgs like events around the 
locations of M H ~ 126 GeV (ATLAS) and M H ~ 124 GeV (CMS) jSHZj. Further searching 
for Higgs boson and studying its properties are still the important tasks for the present and 
upcoming high energy colliders. 

Despite the tremendous success of the SM in describing the high energy physics at the 
energy scale up to several hundred GeV, the instability of the Higgs boson mass leads to 
the "hierarchy problem" [8] which comes from the quadratic loop corrections to the Higgs 
boson mass. In order to give a proper electroweak symmetry breaking (EWSB) scale, the 
Higgs boson mass needs unnatural fine-tuning when it gets a radiative correction with the 
cutoff scale about 10 TeV. In order to solve the "hierarchy problem", physicists developed 
several new particle models such as supersymmetry extra dimensions [10], little Higgs 
[TT| [12] , technicolor and so on. Among these theories, the little Higgs models are proposed 
as one kind of models without fine-tuning in which the Higgs boson is naturally light as a 
result of nonlinearly realized symmetry. The littlest Higgs model (LHM) [T3l [HI [15] is the 
most economical model of them and a phenomenological viable model. 

There are an SU(5) global symmetry and a locally gauged subgroup Gi®G 2 = [SU(2) 1 <S> 
U(1) 1 ]®[SU(2) 2 ®U(1) 2 ] in the LHM. At the scale A s , the global symmetry SU(5) is broken 
into its subgroup 5*0(5). At the same time, the local gauge symmetry [SU(2) £g> U(l)} 2 is 
also spontaneously broken into its diagonal subgroup SU{2) L <g> U(1)y, which is identified as 
the SM gauge group. In the LHM, a set of new heavy gauge bosons (W^, Z H and A H ) and a 
new heavy-vector-like quark (T) are introduced to cancel the quadratic divergence induced 
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by SM gauge boson loops and the top quark loop, respectively. These new gauge bosons 
might provide the significant signatures at the present and future high energy colliders. The 
pp — > Z°H° + X process is one of the main production mechanisms of Higgs boson with 
moderate mass at the LHC, which gives a very distinctive signature. This process could be 
used to measure the Higgs mass and the couplings between Higgs boson and gauge bosons 
and determine the quantum numbers of the Higgs boson. Therefore, investigating the process 
pp -> Z°H° + X at the LHC in the context of the LHM is necessary for probing the LHM 
physics [16]. We find that the SM and minimal supersymmetric standard model (MSSM) 
analyses to the pp — > Z°H° + X process at the LHC have been already existed in Ref . [T7] . 

In this work we study the effects of the LHM on neutral Higgs boson production associated 
with Z° boson up to the QCD next-to-leading order (NLO) at the CERN LHC. In the LHM 
the new neutral gauge bosons, such as Z H and A H , give additional contributions to this 
process. The paper is constructed as follows: In section II, we provide related theory of the 
LHM to our calculations. In section III, we describe the calculations at the leading order 
(LO) and the QCD NLO for the pp — > Z°H° + X process. The numerical results and 
discussions are presented in section IV. Finally, a short summary is given. 



The LHM is based on an SU(5)/SO(5) nonlinear a model. The vacuum expectation value 
(VEV) breaks the SU(5) global symmetry into its subgroup SO (5) and breaks the local 
gauge symmetry [SU(2) ® U{1)} 2 into its diagonal subgroup SU{2) L <g> U{1)y at the same 
time, which is identified as the SM electroweak gauge group. The gauge fields and B ,fl 
associated with the broken gauge symmetries are related to the SM gauge fields by 



II. Related theory of LHM 



W'" = -cW? + sW%, 



W 11 = sW^ + cW 2 M , 



(2.1) 



dB{ + s'B$, 



s'Bf + c'B%, 



(2.2) 
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with mixing angles of 



9i 



9'x 



(2.3) 



At the scale / the SM gauge bosons remain massless, while the heavy gauge bosons acquire 
masses of order /. The W and B are identified as the SM gauge bosons, with couplings of 
g = g\S = g2C and g' = g[s' = g' 2 c'. The electroweak symmetry breaking (EWSB) gives the 
masses for the SM gauge bosons and induces further mixing between the light and heavy 
gauge bosons. We denote the light gauge boson mass eig enstates as W ± (W±), Z°(Z° L ) and 
7(Al) an d the new heavy gauge boson mass eigenstates as W^, Z H and A H . The masses of 
the charged and neutral gauge bosons to the order of v 2 / p are given by 



V 2 f 1 1/2 2^2\ 



(2.4) 
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M 2 = 0, M 2 = m 2 z S 2 w 
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with 



X 



Afv' 



Xh 



BSwCw scs'c'(c 2 s' 2 + s 2 d 2 ) 



(2.9) 



where m z = gv/(2Cw), Cw = cos 6w = ®w is the Weinberg angle, v' and u are the 
VEV's of the scalar SU(2)l triplet and doublet, respectively. In the following numerical 
calculations we take v = 246 GeV and x = 0-5- 
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The couplings of the neutral gauge bosons to quarks are expressed in the form as 
^(1(9 lPl + 9b,Pr) where Pl,r = 1(1 T 7s)- The explicit expressions are given below. 
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9 z r 



9 Z n DD 
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where {7 and D represent the up-type (U = u, c, t) and down-type (D = d, s, b) quarks, 
respectively. The couplings between neutral gauge boson and Higgs boson are expressed as 



9 



hzz = i^vg^ f _ v*_ 
2SICI I f 2 



9 



HZ An — 
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(2.17) 



(2.18) 



The heavy neutral gauge boson Vh (Vh = Zh, Ah) can decay into a fermion pair and 
Z°H°. We obtain the partial decay rates expressed below [T8] . 



12tt L 
T(V h -> 



(^ ff/7 ) 2 (l + 2r f ) + (^ /_/ ) 2 (l - 4r,) ^W f M, 



l v H , 



(9 



V H \2 



1927T 



.V\[(l + r z -r H ) 2 + 8r z ] M V[1} 



(2.19) 
(2.20) 
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where N c = 3 is the color factor, gpff = (g^ff + ^^~)/2, g^ff = (g^ff _ g ]^)/2, 
g A H = g '( c ' 2 - s ' 2 )/(2c's'), g z » = g[c 2 - s 2 )/(2cs), X = l+r 2 +r 2 H -2r z -2r H -2r z r H , and 
Ti = X 2 /My H (Xi — rrif, M z , Mff). Since in our investigated parameter space the Vh —> TT 
and Vh — > Tt(Tt) decays are kinematically forbidden, we assume that the total decay width 
IVff (Vh = Zh, Ah) is the sum of T(Vh — > ff) and T(Vh — > Z°H°), where f = u,d, c, s, b, t, 

III. Analytical calculations 
IIL.l LO calculations 

We generate the Feynman diagrams and their corresponding amplitudes by using Fey- 
nArts3.5 package [19], and apply FormCalc5.4 package [20] to implemented the amplitude 
simplification. The LO contribution to the pp — » Z°H° + X process comes from q-q anni- 
hilation (q = u, d, c, s, b). We denote the partonic process as 

q(pi) + q(P2) -+ Z°(p 3 ) + H°(p 4 ), (q = it, d, c, s, b), (3.1) 

where Pi, P2, P3 and p^ represent the four-momenta of incoming partons, the outgoing Z°- and 
if°-boson, respectively. We use the 't Hooft-Feynman gauge throughout our calculations. 
Comparing with the partonic colliding energy at the LHC, the quark masses, m q (q = 
u,d,c,s,b), are relatively small. We neglect their masses in our further calculations. The 
Feynman diagram for the qq — » Z°H° subprocess in the SM is shown in Figfjja). The 
amplitudes corresponding to FigsJT^a), (b) and (c) without introducing the decay widths in 
propagators are denoted as M z LO (T z = 0), M z l %(T Zh = 0) and Mi%(T AH = 0), respectively. 

As shown in Eq. (I2.10[) -Eq. ( 12TT3]) . the coupling between Z° and quarks in the LHM can 
be obtained from the SM one with a correction of 0(v 2 /f 2 ). The qq — > Z°H° subprocess 
in the LHM obtains additional contributions coming from the diagrams with the exchange 
of heavy gauge bosons Zh and Ah shown in FigJT](b)-(c). These two heavy neutral gauge 
bosons, Zh and Ah, are potentially resonant. For disposal of the singularities due to Vh 
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(a) (b) (c) 

Figure 1: The LO Feynman diagrams for the qq — > Z°H° (q = u, d, s, c, b) partonic process. 

(Vh = Zh, Ah) resonances in the calculations, we have to introduce the decay widths of Zh 
and Ah by doing the following replacements in the resonance propagators of the amplitudes 
M Z L U T z H = 0) for FigfUb) and M^(T Ah = 0) for Figfflc), 

' 1 i75 1 777 , (3-2) 



s 12 - M* h sia - My H + %M Vh T Vh 

where Tv H (Vh = Zh, Ah) represents the decay width of Vh- Then we get the LO amplitudes 
for FigJU(b) and (c) at the tree-level respectively expressed as 

y Sl2 — Ml y 

S12 — M Zh + ijv[ Zh i Zh 
M% = Zl'+h r <o$A H = V). (3.3) 

S\2 — M Ah + IJVI Ah L a h 

The modified amplitudes M. L q and M. L q are safe amplitudes being free of the Zh and Ah 
resonance singularities. Since the 0(a s ) corrections do not contribute to the LO Z H and A H 
decay widths, these replacements cannot induce the double-counting problem in our NLO 
calculations. 

The LO cross section for the subprocess qq — > Z°H° is expressed as 

1 1 (2tt) a C color 

&&> = 4 / E \M L o\ 2 dn 2 , (q = u, d, c, a, b) (3.4) 

spin 

where the factors | and | come from the averaging over the spins and colors of the initial par- 
tons respectively, s is the partonic center-of-mass energy squared, and M.lo is the amplitude 
of all the LO diagrams shown in FigJTJ The summation is taken over the spins and colors of 
all the relevant particles in the qq — > Z°H° subprocess. The integration is performed over 



the two-body phase space of the final particles Z° and H°. dVt 2 is the two-body phase space 
element expressed as 



Within the framework of the QCD factorization, the LO cross section for the process 
pp — > Z°H° + X at the LHC can be obtained by performing the following integration of the 
cross section for the subprocess qq — > Z°H° over the partonic luminosities (see Eq. fl3.6p ). 



where Gi/A(x,fij) {i = u,d,s,c,b) is parton distribution function (PDF) of proton A (= 
Pi, P 2 ) which describes the probability to find a parton i with momentum xpA in proton A, 
s is defined as the total colliding energy squared in proton-proton collision, s = Xix 2 s, and 
lif is the factorization scale. In our LO calculations, we adopt the CTEQ6L1 [2T] PDFs. 

III. .2 Virtual and real emission corrections 

The QCD one-loop vertex correction diagrams for the partonic process qq — > Z°H° with 
nonzero contribution are presented in FigfSJ There exist both ultraviolate (UV) and soft /collinear 
infrared (IR) singularities in the one-loop diagrams. We regularize all the singularities by 
using the dimensional regularization method in D = 4 — 2e dimensions, and apply the modi- 
fied minimal subtraction (MS) scheme to renormalize the relevant fields. The UV divergence 
of the virtual corrections are removed by renormalized wave functions of the relevant quarks. 
We define the renormalization constants of the relevant quark fields as 




(3.5) 




ss,cc,bb, i i 



&LO = 



dXi I dx 2 [G i/Pl (xi, fJ l f)G j /p 2 (x 2 , fJLf) + (Xi <B- X 2 ,P\ <r> Pi)] ctlo(s = XiX 2 s) 

Jo 

(3.6) 




(3.7) 





(1) (2) (3) 

Figure 2: The QCD one-loop vertex correction diagrams for the subprocess qq — > Z°H° (qq = 
uu, dd, cc, ss, bb). 

The notations used in above equation are defined as Ajjv — ~ 1e + ln(47r) and 

A IR = l/ej B - j E + ln(47r). 

Although the total NLO QCD amplitude of subprocess qq — > Z°H° is UV finite after 
performing renormalization procedure, it still contains soft/collinear IR singularities. The 
soft IR singularity can be completely canceled by the contribution of real gluon emission 
subprocess qq — > Z°H°g , while the collinear singularity is eliminated partially by the light- 
quark emission subprocesses q(q)g — > Z°H°q(q) . The remaining collinear IR divergence 
can be absorbed by the counterterms of PDFs. We adopt the analytical expressions for 
IR-singular parts of loop integrals from Ref.[22], and use the expressions in Refs. [23 | |24" | [25] 
to implement the numerical evaluations of IR-safe iV-point(iV < 4) integrals. The Feynman 
diagrams for real gluon/light-quark eimission are depicted in Figj3]and FigJlJ respectively. 

We apply the two cutoff phase space slicing (TCPSS) method [26J to isolate the soft 
and collinear IR singularities of the real emission correction from the IR-safe region. In 
performing the calculations with the TCPSS method, we should introduce arbitrary small 
soft cutoff 6 S and collinear cutoff 6 C . The phase space of the q{pi)q{p2) Z Q {p^)H Q {p4)g{p 5 ) 
partonic process can be split into two regions, E 5 < 5 s \/I/2 (soft gluon region) and E 5 > 
8 s \fi/2 (hard gluon region) by soft cutoff 8 S . The hard gluon region is separated as hard 
collinear (HC) and hard non-collinear (HC) regions by cutoff 8 C . The HC region is the phase 
space where — £15 (or — £25) < 8 c s (ti 5 = (pi — p 5 ) 2 and t 2 5 = (P2 — Ps) 2 )- The phase space of 
light-quark emission q(pi)[(q(pi)]g (P2) — > Z° (p 3 ) H° {pi)q{p^)[q{p^)\ is split into hard collinear 
(HC) region and hard non-collinear (HC) region by introducing a cutoff 5 C . The real gluon 
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Figure 3: The Feynman diagrams for the real gluon emission subprocess qq — > Z°H°g (q 
u, d, c, s, b). 



emission corrections over the HC region are finite and can be calculated numerically with 
general Monte Carlo method [27]. Finally, the cross section for the real emission partonic 
process can be written as 



HC- 



(3.9) 



With the NLO correction components given above, the full QCD NLO corrected cross 
section for the Z°H° production at the LHC can be formally obtained by the QCD factor- 
ization formula as 

a NLO (pp^Z H° + X) = 

j dx Pl dxp 2 |^ [G i/Pl (x Pl , [i^Gj/p^xp,, [if)a% LO {x Pl xp 2 s, fi r )] + (Pi P 2 )| , 

(3.10) 

where i and j run over all possible initial partons contributing to the pp — > Z°H° + X 
process up to the QCD NLO, and the notations of x Pl , x P2 are the same as in Eq. fl3.6p . 
We adopt the CTEQ6m PDFs for Gi/ Pl (x Pl , /i/) and Gj/p 2 (x P2 , p,f) in the NLO calculations 
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(4) (5) (6) 

Figure 4: The Feynman diagrams for the real light-quark emission subprocess q(q)g — > 
Z°H°q(q) (q = u,d,c,s,b). 

[2T] . The total QCD NLO corrected cross section for partonic process qq — > Z°H° can be 
expressed as 

°NLO = °LO + A °NLO =VlO + °R+ & V> ( 3 - 1 1 ) 

where <5"^ ,<5"^ and dy denote the cross sections for tree level, real emission and virtual 
corrections for parton level process, respectively. 

IV. Numerical results and discussions 

In this section we provide and discuss the numerical results for the pp — » Z°H° + X process 
in the LHM up to the QCD NLO. In order to make a cross check with previous work on the 
Z°H° associated production in the SM at the LHC, we take \i = \ij = \i r = y/szn, and the 
input parameters and PDFs being the same as used in Ref.[17j. and calculate the LO and 
NLO QCD corrected total cross sections for pp -»■ Z°H° + X at the y/s = 14 TeV LHC 
in the SM. We get the total cross sections for M H = 140 GeV as a^o = 0.46827(3) pb and 
&nlo — 0.5770(4) pb, separately. The corresponding results can be read out from Table 8 
of Ref . [T7] : (T^o — 0.4684(2) pb and (Jnlo — 0.5768(2) pb, which are coincident with ours 
within the calculation errors. 
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In our following numerical calculations we take the colliding energy in proton-proton 
center-of-mass system as \/s = 8 TeV for the early LHC and y/s = 14 TeV for the future 
LHC. We use one- and two- loop running a s (/x) by taking = 165 MeV and A^ IS = 
226 MeV for the LO and NLO calculations, respectively [2H]- The factorization and the 
renormalization scales are set to be equal for simplicity (/i = [if — u r ). We take // = //q = 
(M H + Mz)/2 in default unless otherwise stated. We neglect the masses of u-, d-, c-, s-, and 
6-quark, and take 

G M = 1.16637 x 1(T 5 GeV- 2 , M w = 80.399 GeV, 

M z = 91.1876 GeV, m t = 171.2 GeV, M H = 125 GeV. (4.1) 

The G M scheme is adopted, i.e., the electromagnetic coupling a is derived from olq = 
\f2G VL My / (1 — Mpp/M§) /it. Considering the constraints of the electroweak precision data 
on LHM parameters [29], we assumed that 0.1 < c < 0.5, 0.1 < c' < 0.9 and 3 TeV < 
f < 7 TeV, and take the representative input LHM parameter set as c = 0.5, d = 0.22 and 
/ = 4 TeV in our numerical calculations if there is no other statement. From Eqs. fl2.5p . 
(12. 6p and (12. 8p with this input parameter set the masses of the heavy gauge bosons Ma h , 
M w ± and M Zh are obtained as 1.461 TeV , 3.025 TeV and 3.025 TeV respectively, where 
the mass values of heavy vector gauge bosons Z\ and W H are beyond the corresponding 
experimental lower mass limits [50] , 

In order to verify the independence of the total NLO QCD corrections on the introduced 
arbitrary cutoff values of 5 S (<5 C ), we depict the Aa^Lo f° r the pp — » uu — > Z°H° + X process 
in the LHM at the a/s = 14 TeV LHC as the functions of 5 S in Figsj5]^a,b), where we take 
c = 0.2, c' = 0.7, / = 2 TeV, 5 C = 5 S /100 and fi = p,Q. The amplified curve for the 
total NLO QCD correction (Aunlo) for the process pp — >■ uu — >■ Z°H° + X is shown in 
FigJ^b). We can see in FigsJ^a,b) that the total QCD correction to the pp — > Z°H° + X 
process does not depend on the arbitrarily chosen value of the cutoffs 5 S and 5 C . The two- 
body correction (Aa^) and three-body correction (Aa^) and the total QCD correction 
(Actnlo = Aa^ + Aa^) for the pp — >■ uu — >■ Z°H° + X process at the LHC are depicted as 
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the functions of the soft cutoff S s in Figsj5]^a). The curve for Ao-jvlo is presented in Figj5]^b) 
together with calculation errors. We adopt also the dipole subtraction (DPS) method to deal 
with the IR singularities for further verification. The A<jnlo results from the DPS method 
including ±lcr statistic errors are plotted as the shadowing region in Figfj^b). It shows that 
the results by using both the TCPSS method and the DPS method are in good agreement. 
In further numerical calculations we adopt the TCPSS method and fix 5 S = 1 x 10 -5 and 
cL = 1 x 1(T 7 . 
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Figure 5: (a) The NLO QCD corrections to the pp — > uu — > Z°H° + X process in the LHM 
at the a/s = 14 TeV LHC as the functions of the soft cutoff 5 S , where we take c = 0.2, 
d = 0.7 and / = 2 TeV, 5 C = 5 S /100 and fi = fiQ. (b) The amplified curve for the NLO 
QCD correction to the cross section Acnlo- 

We show the integrated LO, NLO QCD corrected cross sections and the corresponding 
K-factor {K{fi) = (y NLO {fi)/a LO {fi)) at the = 14 TeV (y/s = 8 TeV) LHC for the 
process pp — > Z°H° + X as the functions of the factorization/renormalization scale (yu/yUo) 
in FigsJHta) (Figsj6|c)), where we set \i = \i T = fif, fi = (M H + M z )/2, c = 0.5, d = 0.22 
and / = 4 TeV. If we define the scale uncertainty for the pp — > Z°H° + X process as 
T) = \ a (v-5vo)-^(v-0-2vo)\ ^ f rom kh e curves in Figsj6]^a,c) we can figure out the corresponding 
uncertainties at the = 14 TeV LHC being rfi% = 0.251, 7] S N % = 0.034, rj^ M = 0.176 and 
V LHM = 0.045, and at the = 8 TeV LHC r?£^ = 0.080, r] s N % = 0.081, r]£g M = 0.033 and 
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rjjtflo = 0.081, respectively. We can see that at the y/s = 14 TeV LHC the LO cross sections 
are strongly related to the scale in the plotted fx range, and the NLO QCD corrections 
significantly reduce the scale uncertainties. But at the y/s = 8 TeV LHC there is no 
distinct improvement for the scale dependence when the QCD NLO corrections are involved. 

Figj6]^b) and Figj6]^d) present the relative deviations defined as 5 (fx) = a y M ^ , as 

the functions of fi/fio, which correspond to Figj6]Ja) and Figj6]Jc), respectively. The two 
figures demonstrate that the NLO QCD corrections obviously reduce the relative deviation 5 
in our plotted fi/fio range. The theoretical relative deviations including the NLO corrections 
are above 9.75% and 2.67% at the y/s = 14 TeV and y/s = 8 TeV LHC, respectively. We can 
read out from the figures that the effects (5) from the heavy neutral gauge boson interactions 
at the y/s = 14 TeV LHC in the vicinity of fx = fi can be about 12.83% for <j LO and 10.37% 
for o nlo . 

In following analysis we show the influence of the LHM parameters c, d , and the global 
symmetry breaking scale /. In FigsJ7(a,b,c,d) we assume fx = fx , c = 0.5 and d = 0.22, and 
depict the plots for the LO and NLO QCD corrected cross sections and the corresponding 
K-factors for the pp — > Z°H° + X process in both the SM and the LHM as the functions 
of the global symmetry breaking scale / at the y/s = 14 TeV and y/s = 8 TeV LHC in 
FigsJTla) and (c), separately. The corresponding relative deviations of the cross sections in 

the LHM from those in the SM, 5(f) = rj^ , are shown in FigsJT^b) and (d), 

respectively. We can see from FigsJ7(a,b,c,d) that when / — > oo, the relative deviations at 
both the LO and the NLO tend to be vanished, and the relative deviations become to be less 
than 5% for the = 14 TeV LHC and the y/s~ = 8 TeV LHC in the ranges of / > 5 TeV 
and / > 4 TeV, respectively. We find also that the deviations are sensitive to the scale / in 
the range of / < 5 TeV in both FigsJT^b) and (d). 

From Eq.f l2.8j) we can conclude that the mass of the heavy gauge boson Zh is mostly 
related with the scale / and the mixing angle parameter c between two SU(2) gauge bosons, 
but not sensitive to the parameter d. For demonstrating the effects from the interactions 
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Figure 6: In these four plots we take c = 0.5, d = 0.22 and / = 4 TeV. (a) The dependence of 
the LO and the QCD corrected cross sections and the corresponding K-factor for the process 
pp — > Z°H° + X on the factorization/renormalization scale (/i//io) at the y/s = 14 TeV^ 
LHC. (b) The corresponding relative deviation of the integrated cross sections in the LHM 
from those in the SM, as the functions of \xj //o at the ^fs = 14 TeV LHC. (c) The LO and the 
QCD corrected cross sections and the corresponding K-factor for the process pp — > Z°H°+X 
versus the scale fi/fio at the y/s = 8 TeV LHC. (d) The corresponding relative deviation of 
the integrated cross sections in the LHM from those in the SM, as the functions of /i//i at 
the v/i = 8 TeV LHC. 
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Figure 7: We take c = 0.5, d = 0.22 and \i = ji^. (a) The LO and NLO QCD corrected cross 
sections and the corresponding K-factors for the pp —> Z°H°+X process at the y/s = 14 TeV 
LHC in both the SM and LHM as the functions of scale /. (b) The relative deviations of the 
cross sections in the LHM from those in the SM corresponding to FigJT^a) as the functions of 
scale /. (c) The LO and NLO QCD corrected cross sections at the \fs = 8 TeV LHC in both 
the SM and the LHM as the functions of scale /. (d) The relative deviations corresponding 
to FigJ7|c) as the functions of scale /. 
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involving Z H boson, we take /i = /i , / = 4 TeV and d = l/y/2 in FigsJH]^a,b,c,d), in which 
case the contributions from the Ah exchange diagrams are vanished (see Eqs. fl2.18p ). We 
plot the LO and NLO QCD corrected cross sections and the corresponding K-factors as the 
functions of the parameter c at the yfs = 14 TeV and \/s = 8 TeV LHC in FigsJHJ^a) and 
(c), separately. The relative deviations of the cross sections in the LHM from those in the 

SM, 5(c) = ffSMM ' are shown in FigsJHJ^b) and (d). FigsJEfa) and (c) show that the 

K-factors in the LHM and SM are beyond 1.29 for both the \/s = 14 TeV and y/s = 8 TeV 
LHC. We can see from Figs]Hb,d) that the difference between the relative deviations of 
Slo(c) and 5nlo(c) goes up with the increment of the mixing angle parameter c in the range 
of c G [0.1, 0.5], and the LO and NLO deviations in the LHM and SM are all sensitive to 
the mixing angle parameter c. We see also that in the range of c < 0.2 the LO and NLO 
relative deviations between the two models are nearly the same for both the \/s = 14 TeV 
and y/s = 8 TeV LHC. 

Eq. (12. 6 1) tells us that the heavy photon mass Ma h mainly depends on the scale / and the 
mixing angle parameter d between two U(l) gauge fields, but is insensitive to the mixing 
parameter c. In order to investigate and discuss the contributions of Ah exchange diagrams 
to the Z°H° associated production, we present the LO and NLO QCD corrected cross 
sections and the corresponding K-factors as the functions of the mixing angle parameter d 
at the \fs = 14 TeV and y/s = 8 TeV LHC in FigsJHJ^a) and (c), separately. In Figsj9]^a,b,c,d) 
we take /x = /Jo, / = 4 TeV and c = l/y/2, in this case there is no contribution from the 
Zh exchange diagrams (see Eq. fl2.18l) ). The corresponding relative deviations of the cross 

\v LHM (c')~a SM (c')\ 

sections in the LHM from those in the SM, 5(d) = a SM (c') > are demonstrated in 

FigsHHb) and (d), respectively. We can see from Figsj9]^a) and (c) that the LO and NLO 
QCD corrected total cross sections in the LHM at the early and future LHC are obviously 
related to the mixing angle parameter d in the range of d G [0.10, 0.65], and the K-factors 
in the LHM are sensitive to d in the range of d G [0.10, 0.65]. Figs|9^b) and (d) demonstrate 
that the difference of #lo( c ') ~ ^vlo( c ') becomes smaller when d increases from 0.40 to 0.65, 
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Figure 8: We take d = l/y/2, f = 4 TeV and // = /i . (a) The LO and NLO QCD 
corrected cross sections and the corresponding K-factors for the pp — > Z°H° + X process 
at the y/s = 14 TeV LHC in both the SM and LHM as the functions of parameter c. (b) 
The relative deviations of the cross sections in the LHM from those in the SM corresponding 
to FigJH] versus parameter c. (c) The LO and NLO QCD corrected cross sections at the 
\fs = 8 TeV LHC in both the SM and LHM as the functions of c. (d) The relative deviations 
corresponding to FigE^c) versus c. 
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while in the range of d G [0.65, 0.90] the NLO and LO relative deviations, Snlo(c') and 
Slo{ c '), have almost the same values for the yfs = 14 TeV and yfs = 8 TeV LHC. 




Figure 9: We take c = l/y/2, f = 4 TeV and y = y . (a) The LO and NLO QCD 
corrected cross sections and the corresponding K-factors for the pp — > Z°H° + X process at 
the a/s = 14 TeV LHC in both the SM and LHM as the functions of the parameter c' . (b) 
The relative deviations of the cross sections in the LHM from those in the SM corresponding 
to FigJHIa) versus parameter d . (c) The LO and NLO QCD corrected cross sections at 
the yfs = 8 TeV LHC in both the SM and LHM as the functions of d. (d) The relative 
deviations corresponding to Figj9]^c) versus parameter d. 

As we know, the final Z° and H° bosons are unstable and can be detected experimentally 
via the subsequential leptonic decays of Z° — > p^ pT and H° — > t + t~ . We employ the 
SM leptonic decay branch ratios of Z Q and H° boson in further numerical calculations, 
i.e., Br(Z° -»■ n + yr) = 3.366% and Br{H° t + t~) = 6.5% [28]. Since the transverse 
momentum distributions of /i + and r + of the process pp — > Z°H° — > y + p~r + r~ + X should 
be the same as those of yT and t~ correspondingly, we present only those of y~ and t~. 
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We depict the LO and QCD NLO corrected transverse momentum distributions of final 
jjT and the corresponding relative deviations of the cross sections in the LHM from those 
in the SM at the y/s = 14 TeV LHC in Figs flW a) and (b) separately, where we take 
M H = 125 GeV, c = 0.5, d = 0.22 and / = 4 TeV. All the curves in FigsfTO^a) go down 
with the increment of the fi~ transverse momentum within the plotted p^ range. The 
differential cross sections, da lo,nlo I dpj. , and the corresponding relative deviations of the 
cross sections at the y/s = 14 TeV LHC in the LHM from those in the SM as the functions 
of p T T are shown in Figs flW c) and (d), respectively. There we adopt again Mh = 125 GeV, 
c = 0.5, c' = 0.22 and / = 4 TeV. In Fig JTUT c) we see that the curves for the LO and 
QCD NLO distributions of p T T in both the SM and the LHM frameworks fall down when 
the transverse momentum p r T goes up. FigsJTUlfa.c) demonstrate that the LO differential 
cross sections of da^ojdp^, and da^o/dp^ in both the SM and the LHM frameworks are 
significantly enhanced by the QCD NLO corrections. FigfTOTb) and FigfTOTd) show that the 
corresponding relative deviations between the two models are significantly suppressed by the 
QCD NLO corrections, and in the ranges of pj, > 130 GeV and p T T > 150 GeV the QCD 
NLO corrected relative deviations can exceed 10%, separately. 

Similar with Figs fiW a.b.c.d) we plot the corresponding distributions of final /i~ and r~ 
in the pp -> Z°H° ->■ ^^t+t- + X process at the y/s = 8 TeV LHC in Figsjma,b,c,d). 
From Fig JTTT a) and Fig JTTT c) we can see that for the y/s = 8 TeV LHC all the curves for 
both the LO and QCD NLO distributions of pj. and p T T decrease with the increment of the 
corresponding transverse momentum, which are similar with the curves for the y/s = 14 TeV 
LHC. Again, we see that both the LO differential cross sections of pt^ and p T T [da^ojdp^ , 
daio/dPr ) are significantly enhanced by the QCD corrections. We can see from Fig JTlT b) 
and Fig JTTT d) that the relative deviations between the two models are significantly suppressed 
by the QCD NLO corrections, and in the ranges of pj, > 155 GeV and p T T > 170 GeV the 
QCD NLO corrected deviations at the y/s = 8 TeV LHC can exceed 10%, separately. 
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Figure 10: The LO and NLO QCD corrected distributions of the transverse momenta of 
final leptons and the corresponding K-factors for the pp — > Z°H° — > {i + fi~r + r~ + X process 
at the a/s = 14 TeV LHC in both the SM and LHM, where we take fi = fio, c = 0.5, 
d = 0.22 and / = 4 TeV. (a) The distributions of pj, . (b) The relative deviations of the 
cross sections in the LHM from those in the SM corresponding to Fig JToT a) as the functions 
of p%, . (c) The distributions of p T T . (d) The relative deviations corresponding to Fig JToT c) 
as the functions of p T T . 
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Figure 11: The LO and NLO QCD corrected distributions of the transverse momenta of 
final leptons and the corresponding K-factors for the pp — > Z°H° — > {i + fi~r + r~ + X process 
at the y/s = 8 TeV LHC in both the SM and LHM, where we take c = 0.5, d = 0.22 and 
/ = 4 TeV. (a) The distributions of pj, . (b) The relative deviations of the differential cross 
sections in the LHM from those in the SM corresponding to Fig fLTT a) as the functions of 
. (c) The distributions of p T T . (d) The relative deviations corresponding to Fig JTlT c) as 
the functions of p T T . 
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V. Summary 

In this paper we investigate the phenomenological effects induced by the new heavy neutral 
gauge bosons in the LHM up to QCD NLO on the Z°H° associated production at the early 
(y/s = 8 TeV) and future (y/s = 14 TeV) LHC. We study the dependences of the LO 
and NLO QCD corrected cross sections on the factorization/renormalization scale fx, the 
LHM parameters c, d and /, and present the LO and NLO QCD corrected distributions of 
the transverse momenta pj, and p T T . It demonstrates that the new neutral gauge bosons 
could induce significant discrepancies to the kinematic observables from the standard model 
predictions for this process at both LO and up to QCD NLO. Our results show that when 
we take the c = 0.5, d = 0.22, / = 4 TeV and fj, = /x , the effects from the heavy neutral 
gauge boson interactions can make the relative deviations to be about 12.83% and 10.37% 
at the LO and up to QCD NLO, respectively. We find that the QCD corrections at the 
y/s = 14 TeV LHC can obviously make the cross section being mildly related to the \x scale, 
and significantly enhance the differential cross sections of the transverse momenta of the 
final decay products \x and r. We also find the LO relative deviations of the integrated cross 
sections are significantly suppressed by the NLO QCD corrections. We conclude that the 
precision measurement of the Z°H° associated production process at the LHC could provide 
the clue of the LHM physics. 
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